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O '. Abstract 
5 ■ 

' A comparative study of the full perturbative QCD calculation of the 

Q^' hadronic production of the Be and B* mesons and the fragmentation approx- 
(D ■ 

' imation is presented. We examined the various subprocesses in detail and the 



distribution C{z), where z is twice the fraction of the Be or B* meson energy 



' in the center of mass of the subprocess, which proved insightful and revealed 

the importance of certain non-fragmentation contributions. We concluded 
that the condition for the applicability of the fragmentation approximation 
in hadronic collisions is that the transverse momentum Px of the produced 
Be and B* be much larger than the mass of the Be meson; i.e., Pt ^ ^Bc- 
Numerical results for the cross sections at the Fermilab Tevatron are also 
presented using updated parton distribution functions with various kinematic 
cuts. 
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The heavy flavored be meson states have attracted considerable interest due to their 
interesting properties. The experimental search for these mesons is now under way at high 
energy colliders such as the LEP e^e~ collider and the Tevatron pp collider. Results are 
expected in the near future. An important theoretical task, relevant to these experiments, is 
to calculate the production cross sections. Since they are weakly bound states hadronization 
is relatively simple and their production cross sections can be accurately predicted in the 
framework of Perturbative QCD (PQCD). At LEP the production of the pseudoscalar ground 
state Be and the vector meson state B* is dominated by decay into a bb pair, followed 
by the fragmentation of a 6 quark into the Be or B* meson Hadronic production, as 

flrst pointed out by Chang and Chen [Q, is dominated at high energies by the subprocess 
gg — s> Bc{B*)bc. The interesting question we address here is how well the fragmentation 
approximation to the full order calculation works for the hadronic production of the Be 
and B* mesons and in which kinematic region. 

The hadronic production of the Be and B* can be calculated fully to order in PQCD, 
which involves a difficult numerical calculation. Since Chang and Chen ^ flrst presented the 
numerical results for the hadronic production, calculations have also been done by several 
other authors p-pHf, not all of which are in agreement. Slabospitsky [Q claimed an order 



of magnitude larger result than Chang and Chen Berezhnoy, et al. |P| obtained a result 
larger than Chang and Chen and smaller than Slabospitsky Masetti and Sartogo |0] 
found a result similar to Berezhony, et al. |^. More recently, Berezhony, et al. ^ found 
that a color factor - was overlooked in their previous work ||5[. After including this factor, 
their revised result is in agreement with Chang and Chen 0. Baranov [^] independently 
also obtained results similar to Chang and Chen 0]. Kolodziej, et al. [|1^ presented results 
using different parton distribution functions and energy scale from that used by Chang and 
Chen so it difficult to directly compare their final result. However, their results for the 
cross sections for the subprocess are similar to others Therefore, we are confident that 
the results of the full order PQCD calculations in Refs. P|^-[I(]| are in agreement. 

An alternative way to calculate the hadronic production is to use the fragmentation ap- 



proximation. From the general factorization theorem it is clear that, for large transverse 
momentum of the or B*, hadronic production is dominated by fragmentation. The 
calculation can then be considerably simplified using this approximation, as was first done 
by Cheung W^. Subsequently, the comparison between the full calculation and the 



fragmentation approximation has been discussed by several authors [0,^|TU|. However, very 
different conclusions have been drawn, although their numerical results are similar. Chang, 
et al. [0 found that when and pr are small the difference between the fragmentation 
approximation and the calculation is large. Berezhony, et al. ^ claimed that the frag- 
mentation approximation breaks down even for very large Pt by examining the ratio of B* 



to Be production. Kolodziej, et al. [|T0| claimed that the fragmentation approximation works 
well if Pt exceeds about 5 — 10 GeV, which is comparable to the Be mass, by investigating 
the Pt distribution of the Be meson. 

It is nontrivial to clarify these points since the full calculation is quite complicated, the 
dominant subprocess involving 36 FejTiman diagrams; but the importance of investigating 
this issue is twofold: From the theoretical perspective, it provides an ideal example to 
quantitatively examine how well the fragmentation approximation works for calculating the 
hadronic production of heavy flavored mesons. In this process both the full contributions 
and the fragmentation approximation can each be calculated reliably. Experimentally, it is 
important to have a better understanding of the production of the B^ mesons in the small Pt 
region where the Be production cross section is the largest. The Pt distribution decreases 
very rapidly as Pt increases. 

We carried out a detailed comparative study of the fragmentation approximation and 
the full ttg QCD perturbative calculation. We first studied the subprocesses carefully. By 
analyzing the singularities appearing in the amplitudes and the Pt distributions of the 
subprocesses we gained very important insight into the processes. We then calculated the 
entire hadronic production cross section for pp (or pp) collisions at the Tevatron energy. To 
facilitate a quantitative comparison we found it very instructive to examine the distribution 
in z, twice the energy fraction carried by the Be meson in the center of mass of the subprocess. 



which is experimentally observable. By investigating the z distribution we found that the 
fragmentation approximation dominates when and only when ^ ^b^- We also found 
that it is insufficient to evaluate the validity of the fragmentation approximation by only 
examining the Pt distribution. In fact, we found that it is simply fortuitous that the Pt 
distribution in the fragmentation approximation for the case of the B^. meson is close to the 
results of the full order a\ calculation when Pt does not satisfy the above condition; i.e., 
Pt ^ Mb^- It is not sufficient for Pt to simply be the same order as M^^, as we discuss 
below. 

In order to obtain some insight into the process we first focus the discussion on the 

subprocess. At the lowest order, cc^, there are 36 Feynman diagrams responsible for the 

dominant gluon fusion subprocess g{k\) + giki) —>■ Bc{p) + b{q2) + c{qi), where ki, k2, p, 

qi, and q2 are the respective momenta. When the energy in the center of mass system, 

V§, is much larger than the heavy quark mass the main contributions to the cross section 

come from the kinematic region where certain of the amplitudes in the matrix element 

are nearly singular; i.e., some of the quark lines or gluon lines are nearly on-shell and the 

related propagators in the Feynman diagrams are nearly singular. This results in the cross 

section for the subprocess coming from the lowest twist contributions being proportional to 
1 /I 

- J' , where Jb^ is the Be decay constant, with some logarithmic correction terms such 



as ln(s/M|^). Now the possible singularities in the square of the matrix element for this 
subprocess must arise from the inverse power(s) of the following factors, or their products, 
which can appear in all the possible denominators of the quark and gluon propagators in 
the Feynman gauge: 



when the Pt of the Be meson is large only (aip + qiY can still be small (~ mf). The 
fragmentation functions can then be extracted from the most singular part containing the 
inverse powers of this factor in the square of the matrix element. It then follows that in the 



qi ■ kj, p ■ kj, (uip + qiY, and {kj - a^p - qiY, 



(1) 



where i,j — 1,2 and ai^2 — 




is the ratio of quark masses. It is easy to see that 
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large Pt region the subprocess is dominated by the fragmentation approximation. However, 
when the Pt of the meson is small the produced Be, as well as the b and the c quarks, can 
be soft or collinear with the beam. In this region the square of the matrix element is highly 
singular because two or more of the internal quarks or gluons in certain Feynman diagrams 
can simultaneously be nearly on-mass- shell. Although this region is a smaller part of the 
phase space, these nearly singular Feynman diagrams, in fact, make a large contribution to 
the cross section and dominate the small Pt region. 

Generally, in the square of the matrix element we can isolate all the terms which con- 



tribute to the lowest twist cross sections using singularity power counting rules |T^. When 
^ Mb^ the lowest twist contributions dominate, while the higher twist contributions 
are suppressed by a factor m^/s and are negligible. We can decompose the terms which 
contribute to the lowest twist cross sections into two components: One component is the 
fragmentation contribution, which dominates the large Pt region. The other is the non- 
fragmentation component, which comes from the other singular parts of the matrix element, 
those in which two or more quarks or gluons are nearly on-shell, as discussed above. This 
non-fragmentation component dominates in the smaller Pt region. The contributions of 
these two components are quite clearly distinguishable in the Pt distribution of the sub- 
process, particularly at large ^/I. In Fig. 1 we show the Pt distribution of the subprocess 
when ^/I = 200 GeV for both the full calculation and the fragmentation approximation 
with as = 0.2, = 4.9 GeV, = 1.5 GeV, and /^^ = 480 MeV. In the fragmentation 
calculation, to reduce the error caused by the phase space integrations, we directly used the 



squared matrix elements [0,0], from which the fragmentation functions are derived [jT],|T4 



rather than the fragmentation functions themselves. It is easy to see in Fig. 1 that when 
Pt is larger than about 30 GeV for the Be and about 40 GeV for the B*, the fragmen- 
tation approximation is close to the full calculation. However, when the value of the 
Pt is smaller than about 30 GeV for the Be and about 40 GeV for the B*, the deviation 
between the fragmentation approximation and the full calculation becomes large and the 
non-fragmentation component clearly dominates their production. This critical value of Pt 



is certainly much larger than the heavy quark masses, or the Be meson mass; we also found 
that it slowly increases with increasing s, which may indicate that there is an additional 
enhancement due to logarithmic terms such as Ins/m^ in the non-fragmentation component 
compared to the fragmentation component. When is not very large this two component 
decomposition is less distinct, since the higher twist terms can not be ignored p|J^,p^. In 
this case, the fragmentation approximation is not a very good approximation, there being 
quite a large discrepancy with the full calculation 0J|,|l3 • 

A similar process is the production of Be and B* mesons in photon-photon collisions, 
which is also instructive to examine more carefully. A comparative study of the full a^a^ 
calculation with the fragmentation approximation in this process was presented in Ref. 
||T^ , p!3[| , where it was claimed that the fragmentation approximation is not valid. There are 
20 Feynman diagrams which can be divided into four gauge invariant subsets corresponding 
to various attachments of photons onto the quark lines; i.e., subsets I, II, III, and IV corre- 
sponding respectively to the attachment of both photons onto the b quark line, onto the c 
quark line, one photon onto the b quark line and the other onto the c, and the interchange of 
b and c. Subset I is dominated by the b quark fragmentation into the Be when Pt is large, as 
discussed above. Subsets III and IV, called recombination diagrams, can only contribute to 
the non-fragmentation component and decrease rapidly when Pt ^ Mb,., as also discussed 
above. However, subset II is somewhat unusual. Although this contribution is relatively 
suppressed by the smaller probability for subsequent bb quark creation, it nevertheless gives 
quite a large contribution to the total cross section because of the enhancement of the c 
quark electric charge. For example, we found that when = 100 GeV the result of the 
full calculation is an order of magnitude larger than the fragmentation calculation, even for 
large Pt- However, when a/S becomes extremely large; e.g., a/J = 800 GeV, the contribu- 
tion of this subset II is dominated by the c quark fragmentating into the Be meson when 
Pt ^ Mb,- This implies that the higher twist terms must have very large coefficients, which 
violates that naive power counting rules which imply that the higher twist terms are sup- 
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pressed by and are therefore negligible when ~ 100 GeV ^ Mb ■ Fortunately, this 

s 

type of contribution, which involves hb quark pair creation, is not important in the hadronic 
production of the Be and B* simply because the factor of 16 enhancement due to the larger 
c quark electric charge is absent. 

Having investigated the subprocess, we next calculated Be and B* production at pp 
colliders, particularly the Fermilab Tevatron. From our above analysis when Pt ^ 
and, therefore, ^ Mb^ the process is h quark fragmentation dominated. However at 
smaller Pt the fragmentation approximation breaks down. In Fig. 2, we compare the Pt 
distributions of the B^. and the B* mesons coming from the full calculation with the frag- 
mentation approximation. In these calculations we used the CTEQ3M parton distribution 



functions |T5|. There is some ambiguity in the choice of the energy scale /i in ^^(/i) and this 
is clearly dependent on the form of the factorization. For example, the component of the h 
fragmentating into the B^ can be factorized with a hh pair first being created at a very short 



distance, with a reasonable choice of this scale being /ii = ym^ + P^^, followed by the h 
fragmentating into a B^ or B*, with a reasonable choice of this energy scale being /i2 = 2mc. 
Factorization in the non-fragmentation component is more complicated and the choice of 
energy scales is not obvious. For simplicity, we will use a uniform choice of energy scales, 
choosing the same scale as set in the h fragmentation component; i.e. ^^(/ii) x a1{^2)- From 
the numerical calculations. Fig. 2, we see that for the B^ meson the Pt distributions for 
Pt > 5 GeV are very similar for the full calculations and the fragmentation approxima- 
tion. This critical value of Pt is much smaller than that found above in the study of the 
subprocess. However, for the B* meson, the result predicted by the fragmentation approxi- 
mation in Fig. 2 differs from that predicted by the full calculation by 50-70% over a much 
larger range of Pt, more consistent with what was found in the study of the subprocess. 
Therefore, it is difficult to decide from only the Pt distributions where the fragmentation 
approximation is reliable for the hadronic production of the B^ and B* mesons. In fact, we 
found that the Pt distribution alone can be misleading. 
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To clarify this issue, we introduce the distribution 

C[z) = j dxidx2g[xi, ii)g{x2, II) ^ , (2) 

where z = (^i ^ and g{xi, //) are the gluon distribution functions. In the subprocess 

s 

center of mass z is simply twice the fraction of the total energy carried by the Be or B* 
meson. The distribution C{z) provides a sensitive means to investigate the dynamics of 
the production process and the fragmentation approximation. Clearly, if the fragmentation 

approximation is valid, " can be factorized as 

dz 

= ^gg^QiQ^ ® Dq^^bX^^ N ■ (3) 

where Dq^^bX^, 1^) are the usual fragmentation functions and ^gg^Q^Q^ is the gluon fusion 
subprocess cross section for production of the heavy quark pair QiQi. In this approximation, 
the integrals over x\ and x^. can be performed, the fragmentation function can be factored 
out, and C{z) is simply proportional to a sum of the usual fragmentation functions which 
is insensitive to the parton distribution functions and to the kinematic cuts. However, if 
the distribution C{z) is quite different from the fragmentation functions, it is an indication 
that the fragmentation approximation is not valid. Therefore, comparing C{z) calculated in 
the fragmentation approximation with the full order calculation, provides a quantitative 
criterion to judge the validity of the fragmentation approximation. We emphasize that z is 
a very useful variable and is an experimentally measurable quantity, at least in principle. 

In Fig 3 we compare the z distribution C{z) calculated in the fragmentation approxi- 
mation with the full order calculation for Be and meson production with a cut of 
Pt > 10 GeV (Fig. 3a) and also with cuts of Pt > 20 GeV and Pt > 30 GeV (Fig. 3b). 
The z distribution C{z) is sensitive to the smallest Pt region for a given P^ cut because the 
Pt distributions of the B^. and mesons decrease very rapidly with increasing Pt- From 
Fig. 3. some general features are evident: For the Be meson, in the fragmentation approxi- 
mation, for smaller Pt cuts the z distributions C{z) in the higher z region is overestimated 
while it is underestimated in the lower z region; but, after integration over z, the result is 



similar to the full calculation. However, for the B* meson, even for the largest Pt cut, 
the z distribution C{z) calculated in the fragmentation approximation is underestimated at 
all values of z and, after integration over z, the result is definitely smaller than the full 
calculation. This feature explains why the Pt distributions shown in Fig. 2 are similar for 
the Be meson even down to Pt ~ Mb^ but are different for the B* meson. This also shows 
that it is simply fortuitous that the Pt distribution of the Be calculated in the fragmentation 
approximation is similar to that from the full calculation for Pt below this value, particu- 
larly down to Pt ~ Mb^- It is also clear that when Pt is increased the distributions become 
closer. As shown in Fig. 3b, when Pt is as large as 30 GeV the curves calculated in the 
fragmentation approximation are quite close to the full calculation. This indicates that 
the fragmentation approximation is valid in the large Pt region, as expected. We emphasize 
here that the difference between the fragmentation approximation and full calculation is not 
universal, but is process-dependent, and thus does not satisfy Atarelli-Parisi evolution. 

Finally, we examined the ratio of cross sections for B* production calculated in the 
fragmentation approximation with the results of the full calculation. As discussed above, 
for B* meson production the fragmentation approximation always underestimates the full 
ttg. The deviation from the full calculation for the B* meson can be used as a criterion 
to test the validity of the fragmentation approximation. The results for the total cross 
section o"(Pt > Pt mm) for various Pt cuts are listed in Table I. Taking agreement within 
30% as the criterion for the validity of the fragmentation approximation we also see that 
Pt should exceed about 30 GeV, a value considerably larger than the heavy quark masses. 
We also compared the ratio of B* to Be production predicted by the full calculation 
with the fragmentation approximation and reached the same conclusion. We note that 
this conclusion is also rather insensitive to the choice the energy scale fi and the parton 
distribution functions. 

Both the spectroscopy and the decays of the B mesons have been widely studied 
|,P!B|-pn|. The excited states below the threshold will decay to the ground state Be 



by emitting the photon(s) or vr mesons. The golden channel to detect the Be meson is 
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Be — > J/ip + T^ip)- However, the branching ratio is quite small, around 0.2. The exclusive 
semileptonic decay mode B^. — > J/ip + 1 + 1/1 has a relatively larger branch ratio, but there 
is "missing energy" . Our calculations should be helpful in searching for the Be since the Pt 
distribution falls so rapidly that the small Pt region is quite important. 

In summary, we carried out a detailed comparative study of the fragmentation approxi- 
mation and the full order QCD perturbative calculation. We first studied the subprocess 
and obtained some insight into the process by analyzing the singularities appearing in the 
amphtude and the Pt distribution of the subprocess. This revealed that there are Feynman 
diagrams present in the full order matrix element in which two or more quarks or gluons 
can be simultaneously nearly on-mass shell and that these can dominate the subprocess over 
the fragmentation approximation in certain regions of the phase space. We then calculated 
the hadronic production process and investigated the z distribution, C{z), which was used 
to test the fragmentation approximation. Prom the study of both the subprocess and the 
hadronic process we conclude that the fragmentation mechanism dominates when and only 
when Pt ^ Mb^ ■ This conclusion is independent of the choice of the energy scales and the 
parton distribution functions. It is only fortuitous that the Pt distribution of the B^. in the 
fragmentation approximation is similar to that of the full calculation for Pt as low as Mb^ ■ 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant DE-PG02-91-ER40684 and by the National Nature Science 
Poundation of China. 
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Table I. Total cross sections o-{Ptb^ > PTmin) in nh for hadronic production of 
the Be and the mesons predicted by the calculation and the fragmentation 
approximation assuming various Pt cuts and |F| < 1.5. The CTEQ3M parton 
distribution functions were used and the values /b^ = 480 MeV, rric — 1.5 GeV, 
rub = 4.9 GeV, and Mb, = 6.4 GeV. 
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(^B*{frag.) 


0.52 


0.55 


0.56 


0.61 


0.63 




0.67 




0.70 
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FIGURES 

Fig. 1. The Pt distributions of the Be and B* meson for the subprocess with — 200 
GeV. The sohd and the doted Hnes correspond to the full calculation and the fragmen- 
tation approximation, respectively. 

Fig. 2. The Pt distributions of Be and 5* meson production at the Tevatron energy 
^/s = 1.8 TeV. The sohd and the doted hnes correspond to the full ctf calculation and the 
fragmentation approximation, respectively. 

Fig. 3. The z distributions C{z) of the B^ and at the Tevatron energy y/s = 1.8 
TeV. The solid lines are the full calculation and the doted lines are the fragmentation 
approximation (a) with the cut Pt > 10 GeV and (b) with the cuts Pt > 20 GeV and 
Pt > 30 GeV. 
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